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Abstract 

This work presents a procedure to extract morphological information from neuronal cells based 
on the variation of shape functionals as the cell geometry undergoes a dilation through a wide 
interval of spatial scales. The targeted shapes are alpha and beta cat retinal ganglion cells which 
are characterized by different ranges of dendritic field diameter. Image functionals are expected 
to act as descriptors of the shape, gathering relevant geometric and topological features of the 
complex cell form. We present a comparative study of classification performance of additive shape 
descriptors, namely Minkowski functionals, and the non-additive multi scale fractal. We found that 
the new measures perform efficiently the task of identifying the two main classes, a and (5, based 
solely on scale invariant information, while also providing intraclass morphological assessment. 
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I. INTRODUCTION 



Many natural phenomena are defined or influenced by the geometrical properties of the 
involved elements, and vice-versa. Examples of such a shape-function relationship include 
the chemical properties of proteins, the aerodynamic efficiency of wings, and the oxygen 
exchanges through elaborated bronchic structures. The close relationship between geometry 
and function provides strong motivation for the geometrical analysis of natural objects. Yet, 
the state-of-the-art of geometrical characterization, an area sometimes called morphometry 
or morphology, remains in a relatively incipient stage where several competing, and often 
divergent, approaches coexist. While powerful methods have been used in physics to express 
relevant geometrical properties, often including differential measurements such as curvature, 
the situation remains particularly challenging in biology. 

The relationship between neuronal shape and function has attracted increasing atten- 
tion due to its far reaching implications for basic neuroscience and for medical applications. 
Neuron morphology has the special characteristic, that it evolves during the developmental 
stage of the cell, being influenced by its molecular environment and the history of synaptic 
activity l|. The mature neuronal shape, together with its membrane electrical properties, 
determine the electric conductance of the cell j3] and account for part of its electrophysiolog- 
ical characteristics such as firing patterns and computational abilities [3,0,0,0]. Software 
packages, such as NEURON, are available for modelling neuronal activity with basis on cable 
theory 0] which can be useful for analysis of real or virtual neurons [s,O,0]- At the same 
time, neuronal shape can vary for different tissues, depending on a number of extracellular 
factors [lillll- Ultimately, it determines the patterns of connectivity and, consequently, the 
overall network computational abilities. On the application side, tasks such as automated 
morphological characterization of neuronal shape and the diagnosis of abnormalities, deserve 
further investigation. 

While the use of morphological tools was severely constrained until recently by the cost of 
the relatively sophisticated systems needed to process images and geometry, the continuing 
advances in computer software and hardware have paved the way for an ever widening 
range of possible applications. Consequently, more effective morphological concepts and 
methods have been developed and reported in the literature, including the use of differential 



geometry concepts such as multi-scale curvature and bending energy 



laliaUJ, methods 
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from mathematical morphology such as skeletonization Q|, as well as the receritlv 
reported framework known as Integral-Geometry Morphological Image Analysis (MIA) [l7|, 



18, 



19[. The latter approach involves the use of additive shape functionals, i.e. mappings 



that take shapes to single scalar values, in terms of a parameter usually related to the spatial 
scale or time. As far as neuroscience is concerned, the contour of a neuronal cell has been 
show, to possess a fractal stn.ctu.e Q a„d ,ts ...tiscale .acta, dimension bee„ used to 
characterize different morphological classes of neuronal cells 2l|. 

Primarily motivated by the possibility of applying the methodology proposed in [l^ 
Q as a novel and potentially useful tool for addressing the problem of neuronal shape 
characterization and classification, the present work also provides an assessment of those 
measures considering a database of real biological data, namely camera lucida images of 
cat ganglion neuronal cells. In order to provide a comparative reference, the multi-scale 
fractal dimension [2^, itself a shape functional, is also considered as a measure for shape 
characterization. 

Integral geometry provides an adequate mathematical framework for morphological image 
analysis, having a core of useful theorems and formulae that in some cases leads to analytical 
results for averages of image functionals 13| while also being quite efficient to implement 
computationally Here, the class of functionals involved are restricted to additive, 

motion invariant and continuous, called Minkowski functionals. These functionals are related 
to usual geometric quantities, for instance, in the Euclidean plane, to area, perimeter and 
connectivity or the Euler number, which expresses the number of holes in a connected pattern 
such as the image of a neuronal cell. 

In order to describe geometrically one object, a set of measures (functionals) is taken and 
the behaviour of these measures is monitored, as some control parameter is varied. In this 
work we compute additive functionals in the plane as the contour of a neuronal cell image is 
inflated by a parallel set dilation of radius r, the control parameter. The non-additive multi 
scale fractal dimension is derived from one of those computed additive functionals, at each 
radius of dilation, giving important [2^ complementary information. 

This article starts by presenting the adopted methodology, the considered shape func- 
tionals, and the statistical procedure for cell identification. The results, which are presented 
subsequently, clearly indicate that the new measures are efficient for distinguishing morpho- 
logically the two functional classes alpha and beta as well as revealing a strong morphological 
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coherence in one of the classes. 



II. METHODOLOGY 

A. Additive shape functionals 

The morphological characterization in Euclidean plane by means of shape functionals 
explores simple properties of convex sets. For these basic geometric objects, such as triangle 
or an ellipse, we may evaluate a change in area while the object undergoes a morphological 
dilation with the knowledge of its initial geometry. For example the change in area of a 
convex body K, after a parallel set dilation using a 2D ball of radius r, can be expressed as 

A{Kr) = A{K) + U{K)r + vrr^, (1) 

where A{K) and U{K) stand for the initial area and perimeter of the object K and r is 
the dilation parameter. The process of taking parallel sets generalises naturally to higher 
dimensions, while the change in hyper volume preserves the general form (P) and is given 
by the Steiner formula 

A^^) = E(!^)w^i'H^K, (2) 



where the coefficients w!)^^ are referred to as quermassintegrals or Minkowski functionals. 



In 

motion invariant and continuous. Moreover, a Theorem by Hadwiger 



7[. These functionals, as a generalisation of known geometric quantities, are additive, 

171 |l9| states that 



these functionals form a complete set of measures, with the above properties, on the set of 
convex bodies 

d 

m) = Y.^^wf{K). (3) 

j=0 

Notwithstanding, the change of any, additive, motion invariant and continuous functional 



can be expressed, using a generalised Steiner formula 
information 

d d-j 
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23[ , in terms of the initial geometric 



</'(^^) = EE-/'/ViSw-'- (4) 

j=0 k=0 ^ ^ 

The notion of connectivity number or Euler characteristic x is central in establishing the 
aforementioned properties of Minkowski functionals. The usual definition of the connectivity 
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from algebraic topology in two dimensions is the difference between the number of connected 
Uc components and the number of holes nh-. 



X{K) =nc- Uh, (5) 

while in three dimensions distinction should be made between two kind of holes namely 
cavities Uhc and handles(tunnels) n^h 

X{K) =nc- nth + Uhc- (6) 

Integral Geometry provides an equivalent definition for connectivity number of a convex set 
K which is given by 

{1 K ^0 
(7) 
K = 0. 

Of great importance is its property of additivity 

xiA) =x(uU^.) = E ^(^*) - E ^(^^ n ir,) + • • • 

i i<j 

... + {-iy+\{K^n...nKi). (8) 



Additivity and motion invariance is inherited by the Minkowski functionals as they are 
related to the connectivity number by the formulae 

lyj^n^) = [ xiAn E,)diiE, z/ = 0, . . . , d - 1 

Jg 

w\'\A) = uMA), oo, = n'/yr{l + d/2). (9) 

In the above expression E^y stands for an //-dimensional plane in R'^. The integral is to be 
taken for all positions, induced by isometrics Q, of E^, weighted by dfi{Ei,), the kinematical 
density which is in turn related to the Haar measure on the group of motions G, see 



23|. 



To sum up, the Minkowski functionals Wu'^\A), as a generalisation of the usual procedure 
for volume determination, counts the number of possible intersection of a //-dimensional 
plane with the domain A. 

If one is to take advantage of the above additivity property all intersection in (jH} must be 
taken into account. When working on a lattice, there is a more expedient route, exploring 
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TABLE I: Minkowski functionals of elementary open bodies which compose a pixel K. 

the discrete nature of the images and the additivity of the Minkowski functionals, which 
consists of a decomposition of the 2D body A into a disjoint collection of interior bodies, 
open edges and vertices. Following the usual nomenclature we denote the interior of a set A 
hy A = A/dA. For an open interior of a n-dimensional body embedded in a (i-dimensional 
Euclidean space there is the following expression for the Minkowski functionals, 

lyW(i) = u = 0,...,d. (10) 

We may then apply additivity and the lack of connectivity of open sets on the lattice to 
determine the functionals for the body whole 

W!;''\V) = Y,Wy\Nm)nmiV), u = 0,...,d. (11) 

m 

Where nm{V) stands for the number of building elements of each type m occurring in the 
pattern V. For a two dimensional space which is our interest for the present neuron images, 
we display in Table Uthe value of Minkowski functionals for the building elements on a square 
lattice of pixels and their direct relation to familiar geometric quantities on the plane. Using 
the information presented in Table H] and equation pij) we have 

A(V)=n2, U{V) = -4n2 + 2ni, xCP) = ^2 - rii + uq. (12) 

So the procedure of calculating Minkowski functionals of a pattern V has been reduced to 
the proper counting of the number of elementary bodies of each type that compose a pixel 
(squares, edges and vertices) involved in the make up of V. 

In section Um we describe typical results for the evaluation of the above presented additive 
functionals using an actual neuron image. The procedure involves the implementation of, 
first an algorithm for the proper parallel set dilation throughout all permitted radius on 
the square lattice and, second, of an algorithm for the calculation of Minkowski functionals 



by counting disjoint building elements based on t 



undertaking the latter is described in detail in 
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le formulae (1121) . An efficient routine for 



B. Multi scale fractal dimension 



As an example of a related non-additive functional, we add to the previous measures the 
multi scale fractal, an approach which has been applied successfully to neuromorphome- 
try 2^ . The notion of multi-scale fractal dimension refers to a quantitative characterization 
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221. In- 



of complexity and the degree of self-similarity at distinct spatial scales, see 
tuitively, the fractal dimension indicates how much the curve extends itself throughout the 
space. As a consequence, more intricate curves will cover the surround space more effectively 
and will display a higher fractal dimension. This quantity is calculated via the derivative of 
the logarithm of the changing interior the neuron cell image undergoes a dilation. 

As such it is immediately derivable from the first additive functional (Area) in a pixel based 
approach as opposed to a curvature approach, see 



C. Implementation 



We have conducted the evaluation of the functionals described above on a 800MHz ordi- 
nary personal computer running Linux. Both the algorithm for exact dilations on the square 



lattice, according to 



functionals 



2J| and the pixel based algorithm for the estimation of the Minkowski 



12| was implemented in SCILAB-2.6. It took approximately 40s to calculate 



the functionals for each dilation radius. 



III. RESULTS 

We start by describing a measure that is perhaps the simplest one but which holds 
important information not only on its own, but through relation with the multi scalar 
fractal dimension described in section Hi Bl In Figure El we show the typical monotonically 
increasing curve of the interior area of a neuron cell as its contour is inflated by a parallel 
set procedure. To capture the gross structure of this measure we calculate the area under 
the curve, sum", and in a size independent manner, the radius at which the area below 
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this curve reaches half of its value at the end of dilation. The fine structure is given by its 
standard deviation, std". 

In Figure m we show a typical curve of the perimeter of the evolving frontier of the neuron 
cell as the contour is inflated by a parallel set procedure. It is important to observe that this 
measure is affected by errors introduced by the discrete nature of the lattice of pixels and 
the low resolution of the original image, which becomes more acute in beta cells because of 
their reduced size. Nonetheless, this effect tends to be less important after the tenth radii 
or so. To capture the gross structure of this measure, we calculate the area under the curve 
sum^ and the size independent radius at which the area below this curve reaches half of 

2 

its overall value. The fine structure is given simply by the standard deviation of the data, 
stdP. 

In contrast to the preceding measure the connectivity or the Euler number of the neuron 
shape as it undergoes the dilation processes is independent of the resolution of the image. It 
is a measure restricted to the topology of the shape counting essentially the number of holes 
at each radius of dilation no matter the holes are perfectly round or not. The alterations 
of connectivity are subtle from step to step, which is reflected by the complex distribution 
of cusps in Figure The vast amplitude of scale for which there is an abrupt change of 
connectivity is a measure of complexity of this type of cell. Figure IHl shows the particular 
behaviour of the connectivity, for samples from the two classes, as the cell shape undergoes 
a dilation. We take the gross information of this measure by extracting the area under the 
interpolated curve, surriinterp, and a monotonicity index given by 

= —4—' (13) 

s + d + p 

where s, d and p counts respectively the number of times that the curve increase, decrease 
or reach a plateau. This index characterizes a perfect monotonically increasing curve when 
its value is 1 and reaching its minimum value for a curve of high variability. This measure 
is designed specially to explore the multi-scale nature of the neuronal complexity. The finer 
structure is captured solely by the standard deviation stddiff of the difference between the 
interpolated curve and the original data. 

As a measure of complexity, the multi-scale fractal dimension has been experimentally 
found to be related to the connectivity of a shape. Although this relationship is not straight- 
forward, there might be a correlation between these measures as commented bellow. For 
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this image functional, we evaluate the maximum fractal dimension, the mean fractality and 
the standard deviation, respectively max, mean and std-^. A typical curve for the multi scale 
dimension is shown in Figure [7| for the same neuron appearing in Figure ^ 

Among all the considered measurements, we found a good separation of alpha and beta 
type cells for the feature space defined by the perimeter half integral radius, R^, and the 

2 

standard deviation for the fine structure of connectivity, stddiff. Figure |H1 shows the obtained 
alpha and beta clustering with both class exhibiting similar dispersion. Another good result 
for morphological characterization was obtained for the connectivity or the Euler number 
of the cell shape. A feature space involving the index of monotonicity and the integral 
of the interpolated connectivity curve is presented in Figure El More efficient than the 
former in separating the two classes, these measures produce a well localised clustering of 
beta cells characterising its geometrical intricacies. This result suggests that alpha class is 
indeed a more homogeneous category while the beta class may have a morphological subclass 
structure. 

Table im shows the correlation coefficients for the twelve measurements considered in this 
work. Of special interest is the correlation between fractality and connectivity, measures that 
have experimentally been found to represent complementary but not redundant measures 
of complexity. Figure ^1 shows a combination of two measures to form a feature space 
which in this case shows a poor separation of classes which is in accordance with the high 
anti-correlation of the involved measures presented in Table IHl Unusually high correlation 
appears between some measures notably as occurring between area and perimeter, suggesting 
a specific tendency that seems to be particular to the type of data (neurons) and not a general 
rule. 



IV. CONCLUSIONS 



The use of additive shape functionals has been recently considered for the characteriza- 
tion of the geometrical properties of several physical objects |l7, Q. The current article 
explored the use of a representative set of such functionals, namely the area, perimeter and 
connectivity, for the characterization of neural shapes represented in terms of a whole set 
of parallel expansions. The multi scale fractal dimension, a non-additive shape functional, 
was also considered clS db standard for comparison. 
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All the adopted shape functionals consist of functions of the dilating radius of each 
parallel body. For the sake of efficience, the following compact subset of global features was 
selected: the area under the functionals, the value of radius where the area reaches its half- 
value, the standard deviation of the functionals, as well as a novel measurements expressing 
the monotonicity of the functionals and the decomposition of the connectivity functional in 
terms of a low and high variation signals. 

All two-by-two combinations of these measures were investigated visually in order to 
identify the combinations of features leading to more pronounced separations between the 
two classes of considered neural cells, namely cat retinal ganglion cells of type alpha and beta. 
The coefficients of correlation of each pair of measures were also estimated and analysed, 
indicating decorrelation between several of the considered features. The obtained results 
confirmed a differentiated potential of each measurement for neural cell clustering, with 
the features derived from the connectivity functional accounting for the best separation 
between classes. However, further investigations and comparison to more established and 
recent methods (such as in j^]) based explicitly on dendritic morphology will be necessary to 
reveal the best realization of the novel methodology proposed in this article. The biological 
implications of such results are that the two type of cells differ in the distribution of holes 
(defined by the respective dendritic arborizations) for different spatial scales. The obtained 
clusters indicated that the alpha cells exhibit less uniform geometrical properties than the 
beta, suggesting the existence of morphological subclasses. 
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FIG. 1: Initial neuron image. Reprinted with permission from B.B. Boycott and H. Wassle. J. 
Physiol(1974) 240 page 402. 



FIG. 2: The same neuron image after the parallel set dilation. 
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FIG. 3: Area (in pixels) as a function of the scale radius (in pixels) for a typical alpha neuronal 
cell given in Figure ^ Subtle information in this graphic is revealed only through the multi-scale 
fractal dimension. 
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FIG. 4: Perimeter (in pixels) as a function of the scale radius (in pixels) for the alpha neuronal 
cell of Figure ^ Note the expected initial decline and a visible fine structure associated with the 
disappearing and less frequently appearing holes. 
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FIG. 5: Connectivity and the difference between data and interpolation for the typical alpha 
neuronal cell, Figure ^ ^ ^ function of parallel set dilation radius (in pixels). 



>. -50 



-100 



10 

Dilation radius 



o • Alpha 
* • Beta 



15 



20 



FIG. 6: Comparison between the connectivity of typical alpha and beta cells showing a different 
range (in pixels) of complexity for the two class of neurons. 
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FIG. 7: Multi scale fractality of a neuronal type alpha cell. This attribute is obtained from the 
area functional and shows the fine structure that is not revealed by that functional alone. Dilation 
radius in pixels. 
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FIG. 8: Clustering of alpha and beta cells based on Ry^ (in pixels) and connectivity measures. 
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FIG. 9: The feature space based on the index of monotonicity and the integral of the connectivity 
(interpolated) curve, for both alpha and beta cells showing a strong clustering of the beta neuron 
class and a much dispersed alpha class. 
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FIG. 10: A not decisive feature space: related measures with high (anti-)correlation. Half integral 
radius in pixels. 
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TABLE II: Correlation coefficients for the set of the measures extracted from Area, Perimeter, Connectivity and Practahty: Standard 
deviations (std"'^'-'',stddiff). Integrals (sum"'*', suminterp), Half radii (R^y^)' monotonicity index (is), mean value (mean) and max value 
(max). Bold face for absolute values of correlation above 0.5. 
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